The Biased Ideal Rectifier 

By W. R. BENNETT 

Methods of solution and specific results are given for the spectrum of the 
response of devices which have sharply defined transitions between conducting 
and non-conducting regions in their characteristics. The input wave consists 
of one or more sinusoidal components and the operating point is adjusted by bias, 
which may either be independently applied or produced by the rectified output 
itself. 

Introduction 

THE concept of an ideal rectifier gives a useful approximation for the 
analysis of many kinds of communication circuits. An ideal rectifier 
conducts in only one direction, and by use of a suitable bias may have the 
critical value of input separating non-conduction from conduction shifted 
to any arbitrary value, as illustrated in Fig. 1. A curve similar to Fig. 1 
might represent for example the current versus voltage relation of a biased 
diode. By superposing appropriate rectifying and linear characteristics 
with different conducting directions and values of bias, we may approximate 
the characteristic of an ideal limiter, Fig. 2, which gives constant response 
when the input voltage falls outside a given range. Such a curve might 
approximate the relationship between flux and magnetizing force in certain 
ferromagnetic materials, or the output current versus signal voltage in a 
negative-feedback amplifier. The abrupt transitions from non-conducting 
to conducting regions shown are not realizable in physical circuits, but the 
actual characteristics obtained in many devices are much sharper than can 
be represented adequately by a small number of terms in a power series 
or in fact by any very simple analytic function expressible in a reasonably 
small number of terms valid for both the non-conducting and conducting 
regions. 

In the typical communication problem the input is a signal which may 
be expressed in terms of one or more sinusoidal components. The output 
of the rectifier consists of modified segments of the original resultant of the 
individual components separated by regions in which the wave is zero or 
constant. We are not so much interested in the actual wave form of these 
chopped-up portions, which would be very easy to compute, as in the fre- 
quency spectrum. The reason for this is that the rectifier or limiter is 
usually followed by a frequency-selective circuit, which delivers a smoothly 
varying function of time. Knowing the spectrum of the chopped input 
to the selective network and the steady-state response as a function of 
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Fig. 1.— Ideal biased linear rectifier characteristic. 
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Fig. 2. — Synthesis of limiter characteristic. 
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frequency of the network, we can calculate the output wave, which is the 
one having most practical importance. The frequency selectivity may in 
many cases be an inherent part of the rectifying or limiting action so that 
discrete separation of the non-linear and linear features may not actually 
be possible, but even then independent treatment of the two processes 
often yields valuable information. 

The formulation of the analytical problem is very simple. The standard 
theory of Fourier series may be used to obtain expressions for the amplitudes 
of the harmonics in the rectifier output in the case of a single applied fre- 
quency, or for the amplitudes of combination tones in the output when two 
or more frequencies are applied. ' These expressions are definite integrals 
involving nothing more complicated than trigonometric functions and the 
functions defining the conducting law of the rectifier. If we were content 
to make calculations from these integrals directly by numerical or mechanical 
methods, the complete solutions could readily be written down for a variety 
of cases covering most communication needs, and straightforward though 
often laborious computations could then be based on these to accumulate 
eventually a sufficient volume of data to make further calculations un- 
necessary. 

Such a procedure however falls short of being satisfactory to those who 
would like to know more about the functions defined by these integrals 
without making extensive numerical calculations. A question of consider- 
able interest is that of determining under what conditions the integrals may 
be evaluated in terms of tabulated functions or in terms of any other func- 
tions about which something is already known. Information of this sort 
would at least save numerical computing and could be a valuable aid in 
studying the more general aspects of the communication system of which 
the rectifier may be only one part. It is the purpose of this paper to present 
some of these relationships that have been worked out over a considerable 
period of time. These results have been found useful in a variety of prob- 
lems, such as distortion and cross-modulation in overloaded amplifiers, 
the performance of modulators and detectors, and effects of saturation in 
magnetic materials. It is hoped that their publication will not only make 
them available to more people, but also stimulate further investigations of 
the functions encountered in biased rectifier problems. 

The general forms of the integrals defining the amplitudes of harmonics 
and side frequencies when one or two frequencies are applied to a biased 
rectifier are written down in Section I. These results are based on the 
standard theory of Fourier series in one or more variables. Some general 
relationships between positive and negative bias, and between limiters and 
biased rectifiers are also set down for further reference. Some discussion is 
given of the modifications necessary when reactive elements are used in the 
circuit. 
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Section II summarizes specific results on the single-frequency biased 
rectifier case. The general expression for the amplitude of the typical 
harmonic is evaluated in terms of a hypergeometric function for the power 
law case with arbitrary exponent. 

Section III takes up the evaluation of the two-frequency modulation 
products. It is found that the integer-power-law case can be expressed in 
finite form in terms of complete elliptic integrals of the first, second, and 
third kind for almost all products. Of these the first two are available in 
tables, directly, and the third can be expressed in terms of incomplete 
integrals of the first and second kinds, of which tables also exist. No direct 
tabulation of the complete elliptic integrals of the third kind encountered 
here is known to the author. They are of the hyperbolic type in contrast 
to the circular ones more usual in dynamical problems. Imaginary values 
of the angle would be required in the recently published table by Heuman . 

A few of the product amplitudes depend on an integral which has not 
been reduced to elliptic form, and which is a transcendental function of two 
variables about which little is known. Graphs calculated by numerical 
integration are included. 

The expressions in terms of elliptic integrals, while finite for any product, 
show a rather disturbing complexity when compared with the original 
integrals from which they are derived. It appears that elliptic functions 
are not the most natural ones in which the solution to our problem can be 
expressed. If we did not have the elliptic tables available, we would prefer 
to define new functions from our integrals directly, and the study of such 
functions might be an interesting and fruitful mathematical exercise. 

Solutions for more than two frequencies are theoretically possible by the 
same methods, although an increase of complexity occurs as the first few 
components are added. When the number of components becomes very 
large, however, limiting conditions may be evaluated which reduce the 
problem to a manageable simplicity again. The case of an infinite number 
of components uniformly spaced along an appropriate frequency range has 
been used successfully as a representation of a noise wave, and the detected 
output from signal and noise inputs thus evaluated 2 . The noise problem 
will not be treated in the present paper. 

I. The General Problem 
Let the biased rectifier characteristic, Fig. 1, be expressed by 

/ 0, E <b\ 

l = ( (1-1) 

\f(E -b), b< E/ 

1 Carl Heuman, Tables of Complete Elliptic Integrals, Jour. Math, and Physics, Vol. 
XX, No. 2, pp. 127-206, April, 1941. 

2 W. R. Bennett, Response of a Linear Rectifier to Signal and Noise, Jour. Acous. Soc. 
A mer., Vol. 15, pp. 164-172, Jan. 1944. • . . , 
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Then if a single frequency wave defined by 

E = P cos pt, - P < b < P, (1.2) 

is applied as input, the output contains only the tips of the wave, as shown 
in Fig. 3. It is convenient to place the restrictions on P and b given in 
Eq. (1.2). The sign of P is taken as positive since a change of phase may 
be introduced merely by shifting the origin of time and is of trivial interest. 
If the bias b were less than — P } the complete wave would fall in the con- 
ducting region and there would be no rectification. If b were greater than 




Fig. 3. — Response of biased rectifier to single-frequency wave. 

P, the output would be completely suppressed. Applying the theory of 
Fourier series to (1.1) and (1.2), we have the results 



/ = -=■ + zl a n cos n pt 

I n-l 



7T Ja 



arc cos b/P 



f(P cos x — b) cos nx dx 



(1.3) 



(1.4) 



When two frequencies are applied, the output may be represented by a 
double Fourier series. The typical coefficient may be found by the method 
explained in an earlier paper by the author 1 . The problem is to obtain the 
double Fourier series expansion in x and y of the function g(x,y) defined by: 

(0, P cos x + Q cos y < b \ 

(1.5) 
f(P cos x + Q cos y - b), b < P cos x + Q cos y) 

3 VV. R. Bennett, New Results in the Calculation of Modulation Products, B.S.T.J., 
Vol. XII, pp. 228-243, April, 1933. 
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We substitute the special values * = pt,y = qt after obtaining the expansion. 
Let 

k, = Q/P, k = -b/P (1.6) 

The mcst general conditions of interest are comprised in the ranges: 

< *i < 1, - 2 < h < 2 (1.7) 
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Fig. 4. — Regions in ry-plane bounded by £o + cos x + ki cos y = 0. 

The regions in the .ry-plane in which g(x,y) does not vanish are bounded 
by the various branches of the curve : 

k + cos x + £i cos y = (1 .8) 

We need to consider only one period rectangle bounded by x = ±v, y = ±ir, 
since the function repeats itself at intervals of 2ir in both x and y. The 
shape of the curve (1.8) within this rectangle may have three forms, which 
are depicted in Fig. 4. In Case I, £ + h >k, k - h < 1, the curve 
divides into four branches which are open at both ends of the x- and y-axes. 
In Case (2), h + *i <1, h - h > -1, the curve has two branches open 
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at the ends of the y-axis. In Case (3), —1 < k -f- k x < 1, k — ki < — 1, 
a single closed curve is obtained. The limits of integration must be chosen 
to fit the proper case. The Fourier series expansion of g(x,y) may be 
written : 

00 00 . 

g(x, y)=EZ dmn cos mx cos ny (1.9) 

m=0 n-0 

where a m „ is found from integrals of the form: 

A = -^r / dy I f(P cos x "r* Q cos y — b) cos mx cos ny dx (1.10) 

Here, as usual, c m is Neumann's discontinuous factor equal to two when m 
is not zero and unity when m is zero. The values of the limits for the dif- 
ferent cases are : 
Case I, a m „ = A\ + A 2 

(xi = 0, Xt = arc cos (— ko — k\ cos y)' 

Ai = A with limits I I - k I (1-1 D 

- cos — , y-i -■ 

Xi = 0, Xi = T 




1 - *. d-12) 
y2 = arc cos — r— 



Xi = arc cos ( — k — ki cos y 
y* = * 

X2 = arc cos (— ko — ki cos y) 
y 2 = arccos^-^-j 



(1.13) 



(1.14) 



For a considerable variety of rectifier functions/, the inner integration may 
be performed at once leaving the final calculation in terms of a single definite 
integral. 

A somewhat different point of view is furnished by evaluating the integral 
(1.4) for the biased single-frequency harmonic amplitude, and then replacing 
the bias by a constant plus a sine wave having the second frequency. When 
each harmonic of the first frequency is in turn expanded in a Fourier series 
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in the second frequency, the two-frequency modulation coefficients are ob- 
tained. Some early calculations carried out graphically in this way are 
the source of the curves plotted in Figs. 18 to 21 inclusive, for which I am 
indebted to Dr. E. Peterson. 

If reactive elements are used in the rectifier circuit, the voltage across the 
rectifying element may depart from the input wave shape applied tp the 
complete network. The solution then loses its explicit nature since the 
rectifier current is expressed in terms of input voltage components which in 
turn depend on voltage drops produced in the remainder of the network 
by the rectifier currents. Practical solutions can be worked out when 
relatively few components are important. 




Fig. 5 — Biased rectifier in series with RC network. 

As an example consider the familiar case of a parallel combination of 
resistance R and capacitance C in series with the biased rectifier, Fig. 5. 
If C has negligible impedance at all frequencies of importance in the rectifier 
circuit except zero, we may assume that the voltage across R is constant and 
equal to I R, where I is the d-c. component of the rectifier current. The 
voltage across the rectifier unit is then E — I Q R. The effect is a change 
in the value of bias from b to b + hR- If the d-c component in the output 
is calculated for bias b + I R, we obtain the value of J in terms of b + I R, 
an implicit equation denning J . If this equation can be solved for 7 , the 
bias b + I R can then be determined and the remaining modulation products 
calculated. 

A more important case is that of the so-called envelope detector, in which 
the impedance of the condenser is very small at all frequencies contained in 
the input signal, but is very large at frequencies comparable with the band 
width of the spectrum of the input signal. These are the usual conditions 
prevailing in the detection of audio or video signals from modulated r-f or 
i-f waves. The solution depends on writing the input signal in the form. 
of a slowly varying positive valued envelope function multiplying a rapidly 
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oscillating cosine function. That is, if the input signal can be repre- 
sented as 

E = A (/) cos (t), (1.15) 

where A (t) is never negative and has a spectrum confined to the frequency 
range in which 2-kJC is negligibly small compared with \/R, while cos <j>{t) 
has a spectrum confined to the frequency range in which \/R is negligibly 
small compared with 2irfC, we divide the components in the detector output 
into two groups, viz.: 

1. A low-frequency group It/ containing all the frequencies comparable 
with those in the spectrum of A (t). The components of this group flow 
through R. 

2. A high-frequency group Ih/ containing all the frequencies comparable 
to and greater than those in the spectrum of cos <f> (/). The components 
of this group flow through C and produce no voltage across R. 

The instantaneous voltage drop across R is therefore equal to Ii/R, and 
hence the bias on the rectifier is b + h/R- If A and <f> were constants, we 
could make use of (1.3) and (1.4) to write: 



2 [ 

IT JQ 



III + hf m -x +2 a n cos nd (1.16) 

arocos [(6+/j/K)/A] 

f(A cos x — b — Ii/R) cos nx dx (1.17) 



If A and <f> are variable, the equation still holds provided I t /R < A at all 
times. Assuming the latter to be true (keeping in mind the necessity of 
checking the assumption when I u is found), we note that terms of the form 
a n cos » 6 consist of high frequencies modulated by low frequencies and hence 
the main portion of their spectra must be in the high-frequency range. 
Hence we must have as a good approximation when the envelope frequencies 
are well separated from the intermediate frequencies, 

**/= T = " / M cosx - b - I lf R)dx (1.18) 

I TT JO 

This equation defines /// as a function of A, and if it is found that the 
condition b + hjR < A is satisfied by the resulting value of I lf , the problem 
is solved. If the condition is not satisfied, a more complicated situation 
exists requiring separate consideration of the regions in which b + Ii/R < A 
and b + I lf R > A . 

To be specific, consider the case of a linear rectifier with forward con- 
ductance a = l/R, and write V = h/R. Then 

^ V = \/A - (b + V)* - (b + V) arc cos b +V. (1.19) 
K A 
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When 6 = (the case of no added bias), this equation may be satisfied by 
setting 

V = cA,Q<c < 1, (1.20) 



which leads to 



tR /l 

r ' ' y c 2 



- — 1 — arc cos c, (1.21) 



defining c as a function of R /R. The value of c approaches unity when 
the ratio of rectifier resistance to load resistance approaches zero and falls 
off to zero as R /R becomes large. The curve may be found plotted else- 
where 4 . This result justifies the designation of this circuit as an envelope 
detector since with the proper choice of circuit parameters the output 
voltage is proportional to the envelope of the input signal. 

The equations have been given here in terms of the actual voltage applied 
to the circuit. The results may also be used when the signal generator 
contains an internal impedance. For example, a nonreactive source inde- 
pendent of frequency may be combined with the rectifying element to give a 
new resultant characteristic. If the source impedance is a constant pure 
resistance r throughout the frequency range of the signal input but is 
negligibly small at the frequencies of other components of appreciable size 
flowing in the detector, we assume the voltage drop in r Q is r ai cos <£ (/). 
We then set n = 1 in (1.17) and replace Oi by (A — A)/r , where A Q is 
the voltage of the source. The value of I if in terms of A from (1.18) is 
then substituted, giving an implicit relation between A and Ao • 

A further noteworthy fact that may be deduced is the relationship be- 
tween the envelope and the linearly rectified output. By straightforward 
Fourier series expansion, the positive lobes of the wave (1.15), may be 
written as: 

/E, £>0\ p 

Er=[ =A{t) - +-cos*(/) 

\0, £<0/ L7r 2 

_ 2 f, (-) w cos 2m 0(0 1 (1 n) 

fS" 4« 2 -l J 
Hence if we represent the low-frequency components of E T by E lf , we have: 

E„ = 42 (1.23) 

or 

A (/) - vEu (1.24) 

* See, for example, the top curve of Fig. 9-25, p. 31 1, H. J. Reich, Theory and Applica- 
tions of Electron Tubes, McGraw-Hill, 1944. 
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Equation (1.23) expresses the fact that we may calculate the signal com- 
ponent in the output of a half-wave linear rectifier by taking 1/tt times the 
envelope. Equation (1.24) shows that we may calculate the response of 
an envelope detector by taking t times the low-frequency part of the 
Fourier series expansion of the linearly rectified input. Thus two procedures 
are in general available for either the envelope detector or linear rectifier 
solution, and in specific cases a saving of labor is possible by a proper choice 
between the two methods. The final result is of course the same, although 
there may be some difficulty in recognizing the equivalence. For example, 
the solution for linear rectification of a two-frequency wave P cos pt + Q 
cos ql was given by the author in 1933 3 , while the solution for the envelope 
was given by Butterworth in 1929 5 . Comparing the two expressions for 
the direct-current component, we have: 

Eis — — [2E — (1 — k 2 ) K], where K and E are complete elliptic integrals 
w 2 

of the first and second kinds with modulus k = Q/P 

A (t) = — (1 + k) Ei, where Ei is a complete elliptic intregal of the 

second kind with modulus £i = 2 y/kf{\ + *)■ Equation (1.24) implies 
the existence of the identity 

(1 + k) Ei = IE - (1 - k 2 ) K (1.25) 

The identity can be demonstrated by making use of Landen's transforma- 
tion in the theory of elliptic integrals. 

2. Single-Frequency Signal 

The expression for the harmonic amplitudes in the output of the rectifier 
can be expressed in a particularly compact form when the conducting part 
of the characteristic can be described by a power law with arbitrary ex- 
ponent. Thus in (1.4) if/(s) — az v , we set X = b/P and get 

• arc cos X 



a„ = / (cos x — \) cos nx dx 

■K JO 

2*r(, + \)aP'(\ -X)" +1 



-t- i)or u - a; , a i _ \\ 

Tlr( W!J 'G- + *!-*' + 5 ;L t- x ) ™ 

6 S. Butterworth, Apparent Demodulation of a Weak Station by a Stronger One 
Experimental Wireless, Vol. 6, pp. 619-621, Nov. 1929. 
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The equation holds for all real values of v greater than —1. The symbol 
F represents the Gaussian hypergeometric function 6 : 

rv u \ . . ab ,a(a + 1) b (b + 1) 2, (1 9 \ 

The derivation of (2.1) requires a rather long succession of substitutions, 
expansions, and rearrangements, which will be omitted here. 

When v is an integer, the hypergeometric function may be expressed in 
finite algebraic form, either by performing the integration directly, or by 
making use of the formulas: 

FQi/2, — it/2; 1/2; z) = cos (jx arc sin z), 

(2.3) 
/l + /* 1 — ft 3 2 \ _ sin (n arc si n z) 

together with recurrence formulas for the F-f unction. When v is an odd 
multiple of one half, the F-function may be expressed in terms of complete 
elliptic integrals of the first and second kind with modulus [(1 — X)/2] by 
means of the relations, 



2s 2 



FihhUr) =-K, 



7T 



F(-hhUk 2 ) =-E, 



(2.4) 



and the recurrence formulas for the /^-function. For the case of zero bias, 
we set X = 0, and applv the formula 

rft) r('4Q 

F(a, I - .; ,i 1/2) = , c + \( - > 1 + t r ^ ( 2 - 5 > 

obtaining the result : 

2 , r(y+1)r (2^±j) r (2z±i^ f , 

We point out that the above results may be applied not only when the 
applied signal is of the form P cos pt with P and p constants, but to signals 

For an account of the properties of the hypergeometric function, see Ch. XIV of 
Whittakcr and Watson, Modern Analysis, Cambridge, 1940. A discussion of elliptic 
integrals is given in Ch. XXII of the same book. 



a„ = 



(2.6) 
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in which P and p are variable, provided that P is always positive. We thus 
can apply the results to detection of an ordinary amplitude-modulated wave 
or to the detection of a frequency-modulated wave after it has passed through 
a slope circuit. 

A case of considerable practical interest is that of an amplitude-modulated 
wave detected by a diode in series with a parallel combination of resistance 
R and capacitance C. The value of C is assumed to be sufficiently large so 
that the voltage across R is equal to the <io/2 component of the current 
through the diode multiplied by the resistance. This is the condition for 
envelope detection mentioned in Part 1. The diode is assumed to follow 
Child's law, which gives v = 3/2. We write 

/o _ _ _ ao/2 (2^r(3) V' h ' ~~2~7 ( } 

where X = V/P. Note that V is a constant equal to the direct-voltage 
output if P is constant. If P varies slowly with time compared with the 
high-frequency term cos pt, V represents the slowly varying component of 
the output and hence is the recovered signal. 
But 

F(h b 3; k 2 ) = J£ 4 [2(2* 2 - 1)E + (2 - 3k 2 )(l - k 2 )K) (2.8) 

where K and E are complete elliptic integrals of the first and second kind 
with modulus k. Hence 



where the modulus of A' and E is \/(l — X)/2- This equation defines p 
as a function of X, and hence by inversion gives X as a function of p. The 
resulting curve of X vs. p is plotted in Fig. 6 and may be designated as the 
function X = g (p). If we substitute X = V/P we then have 

V = P g (3v/Ra V2P) (2.10) 

This enables us to plot V as a function of P, for various values of Ra, Fig. 7. 
Since P may represent the envelope of an amplitude-modulated (or diff- 
erentiated FM) wave, and V the corresponding recovered signal output 
voltage, the curves of Fig. 7 give the complete performance of the circuit 
as an envelope detector. In general the envelope would be of form P = 
Poll + c s(l)], where s(t) is the signal. We may substitute this value of P 
directly in (2.10) provided the absolute value of c s(t) never exceeds unity. 
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Fig. 6.— The Function X = g( P ) denned by Eq. (2.9). 
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Fig. 7. — Performance of 3/2 — power-law rectifier as an envelope detector with low-imped- 
ance signal generator. 

To express the output in terms of a source voltage P in series with an 
impedance equal to the real constant value r at the signal frequency and 
zero at all other frequencies, we write 



Po - P 



= ai = 



3aP 3/2 (l - 



Sf*'(*-*»» 4*) < 2 



ii) 
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or 



where 



*-(i+S*)* 



(2.12) 



H = 



3Ra(l - \ypl 



- \) 2 j>> / 1 - x \ 

V2 ' F [*»-*'' 3m ' 2 ) 



= ~ V2P[2(1 -k 2 + k A )E - (2 - * 2 )(1 - k 2 )K). 

5w 



(2.13) 
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Fig. 8. — Performance of 3/2 — power-law rectifier as an envelope detector with impedance 
of signal generator low except in signal band. 

By combining the curves of Fig. 7 giving V in terms of P with the above 
equations giving the relation between P and P , we obtain the curves of 
Figs. 8, 9, 10, giving 7asa function of P . The curves approach linearity 
as Ra is made large. On the assumption that the curves are actually linear, 
we define the conversion loss D of the detector in db in terms of the ratio 
of maximum power available from the source to the power delivered to the 
load: 



Z, = 1 O,og lo ^=lOlog lo (0 2 ^ 



(2.14) 



Curves of D vs r /R are given in Figs. 11 and 12. The optimum relation 
between r and R when the forward resistance of the rectifier vanishes has 
long been known to be r /R = .5. The curves show a minimum in this 
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region when Ra is large. In the limit as Ra approaches infinity, we may 
show that the relation between P and V approaches: 



ft-r(i + $) 



(2.15) 
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Fig. 9. — Performance of 3/2 — power-law rectifier as an envelope detector with impedance 
of signal generator low except in signal band. 




Fig. 10. — Performance of 3/2 — power-law rectifier as an envelope detector with impedance 
of signal generator low except in signal band. 



The corresponding limiting formula for D is 



(2.16) 
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The minimum value of D is then found to occur at r = R/2 and is zero 
db. We note from the curves that the minimum loss is 1.2 db when Ra = 
10 and 0.4 when Ra = 100. 

This example is intended mainly as illustrative rather than as a complete 
tabulation of possible detector solutions. The methods employed are 
sufficiently general to solve a wide variety of problems, and the specific 
evaluation process included should be sufficiently indicative of the proce- 
dures required. Cases in which various other selective networks are asso- 
ciated with the detector have been treated by Wheeler 7 . 
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Fig. 1 1. — Conversion loss of 3/2 — power-law rectifier as envelope detector with impedance 
of signal generator low except in signal band. 
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Fig. 12. — Conversion loss of 3/2 — power-law rectifier as envelope detector with impedance 
of signal generator low except in signal band. 

3. Two-Frequency Inputs 

The general formula for the coefficients in the two-frequency case depends 
on a double integral as indicated by (1.10). In many cases one integration 
may be performed immediately, thereby reducing the problem to a single 
definite integral which may readily be evaluated by numerical or mechanical 

7 H. A. Wheeler, Design Formulas for Diode Detectors, Proc. I. R. E., Vol. 26, pp. 
745-780, June 1938. 
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means. It appears likely in most cases that the expression of these results 
in terms of a single integral is the most advantageous form for practical 
purposes, since the integrands are relatively simple, while evaluations in 
terms of tabulated functions, where possible, often lead to complicated 
terms. Numerical evaluation of the double integral is also a possible method 
in cases where neither integration can be performed in terms of functions 
suitable for calculation. 

One integration can always be accomplished for the integer power-law 
case, since the function / (P cos x + Q cos y — b) in (1.12) then becomes a 
polynomial in cos x and cos y. Cases of most practical interest are the 
zero-power, linear, and square-law detectors, in which /(z) is proportional 
to z°, z 1 , and z respectively. The zero-power-law rectifier is also called a 
total limiter, since it limits on infinitesimally small amplitudes. We shall 
tabulate here the definite integrals for a few of the more important low-order 




Fig. 13. — Response of biased total limiter to two-frequency wave. 

coefficients. To make the listing uniform with that of our earlier work, we 
express results in terms of the coefficient A mn , which is the amplitude of the 
component of frequency mp ± nq. The coefficient A m „ is half of a mn when 
neither m nor n is zero. When m or n is zero, we take A mn = dmn and drop 
the component with the lower value of the ± sign. When both m and n 
are zero, we use the designation Aoo/2 for aoo, the d-c term. In the tabula- 
tions which follow we have set /(z) = az" with v taking the values of zero 
and unity. 

We first consider the biased zero-power-law rectifier or biased total 
limiter. This is the case in which the current switches from zero to a 
constant value under control of two frequencies and a bias as illustrated 
by Fig. 13. The results are applicable to saturating devices when the 
driving forces swing through a large range compared with the width of the 
linear region. It is also to be noted that the response of a zero-power-law 
rectifier may be regarded as the Fourier series expansion of the conductance 
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of a linear rectifier under control of two carrier frequencies and a bias. 
The results may therefore be applied to general modulator pioblems based 
on the method described by Peterson and Hussey 8 . We may also combine 
the Fourier series with proper multiplying functions to analyze switching 
between any arbitrary forms of characteristics. We give the results for 
positive values of k . The corresponding coefficients for — k can be ob- 
tained from the relations: 



loo -^00 



(3.1) 



2 2 

A~ n = (-) m+n+i AL, m + n>0 
Here we have used plus and minus signs as superscripts to designate co- 
efficients with £ias +k and - k respectively. We thus obtain a reduction 
in the number of different cases to consider, since Case III consists of nega- 
tive bias values only, and these can now be expressed in terms of positive 
bias values falling in Cases I and II. It is convenient to define an angle 
by the relations: 

e /arc cos Lz_*_° , * + *i > 1, *o - h < 1 (Case I) \ ^ 

\ ,• *o + h < 1, k Q - h > -1 (Case 11)/ 

Zero-Power Rectifier or Total-Limiter Coefficients 
Settingy(z) = a in (1.10), 

i^ = 1 - 1 [ ' arc cos (*„ + h cos y) dy 

2a ir Je 

— = 4 f Vl - (*o + h cos y) 2 dy 



An = 2*i r 

a ir 2 Je 



sin 2 y dy 



Vl — (ko + ki cos y) 2 

dl 1 - 1 f cos Vl - % + ki cos yj 2 dy 

a 7T" J e 

— = ~* f (*n + *i cos y) Vl - (*o + *i cos y) 2 dy 
a ir' Je 

A m = 2*i f r sin 2 y cos y dy 

a ir 2 Je 

A 



Vl — (k + k, cos y)- 

i** = -1 ( (k + k x cos y) cos y Vl - (*o + *i cos y) 2 rfy 
a iH Je 



{3.3) 



8 E. Peterson and L. W. Hussey, Equivalent Modulator Circuits, B. S. T. J., Vol. 18, 
pp. 32-48, Jan. 1939. 
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Similarly 


for 


a linear rectifier: 




* 


A M 

2 


= ^-° + ab 






A To 


= aP - Ato 






Aqi 


= aQ - Ati 






J*mn 





m + n > 1 



C3.4J 



We have shown in Fig. 2 how an ideal limiting characteristic, which trans- 
mits linearly between the upper and lower limits, may be synthesized from 
two biased linear rectification characteristics. Equation (3.4) shows how 
to calculate the corresponding modulation coefficients, when the coefficients 
for bias of one sign are known. The limiter characteristic is equal to az— 
h (z) - h (z), where 

/z-bi, z> h\ I 0, 2 > -b\ 

fi(z) = a[ , /i («) = «( ? (3.5) 

\ 0, z < V \z + h, z < -b 2 / 

The expression for / 2 (2) may also be written : 

/z - (-b 2 ), z> -b 2 \ 
/a (2) = a (z + b 2 ) - a ) (3.6) 

\ 0, 2 < -b 2 J 

Hence the modulation coefficient A mn for the limiter may be expressed in 
terms of A m „ (h) and A m „ (—b 2 ) as follows: 

A mn = -A mn (61) + (-) m+n A mn {b 2 ), m + n ^ 1 (3.7) 

If the limiter is symmetrical {b\ = b 2 ), the even -order products vanish and 
the odd orders are doubled. The terms aP, aQ are to be added to the 
dexter of (3.7) for A i ,A 01 respectively. The odd linear-rectifier coefficients, 
when multiplied by two, thus give the modulation products in the output 
of a symmetrical limiter with maximum amplitude k , as may be seen by 
substituting bi = b 2 = —ko in (3.7). For the fundamental components 
aP and aQ respectively must be subtracted from twice the A 10 and A 01 co- 
efficients for & . 

Linear Rectifier Coefficients 
D.C. 

^/aP = k + -f [Vl - (ko + kicosy)* 
2 *" J « (3.8) 

— (&o + ki cos y) arc (cos £0 + £1 cos y)] dy 
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Fundamentals 
Ai /aP = 1 + -. [ [(*» + ki cos y) Vl - (k + h cos y) 2 

IT* JO 

— arc cos (k + &i cos y)] <fy 

^oi/a/ 3 = *, + A ( [ Vl - (*o + ki cos y) 2 

— (h + h cos y) arc cos (k + h cos y)] cos y </y 

Sum and Difference Products — Second Order 

An = —r I [(k + ki cos y) Vl - (k + £1 cos y) 2 

— arc cos (&o + £i cos y)] cos y dy 
Sum and Difference Products — Third Order 
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(3.9) 



(3.10) 



(3.11) 



2aP 



j [1 - (*o + ki cos y) 2 ] 3 ' 2 cos y dy (3.12) 



37T 2 J fl 

The above products are the ones usually of most interest. Others can 
readily be obtained either by direct integration or by use of recurrence 
formulas. The following set of recurrence formulas were originally derived 
by Mr. S. O. Rice for the biased linear rectifier: 

In A mn + ki (m — m — 3) A m+hn -i 

+ kl (m + » + 3)/l TO+ i,«_l + 2£ » ^m+l.n = 

2» A mn + fa (» + m — 3) i4 m _i,„+i 

+ ^i (« — w + 3)^4 m _i, n+i + 2k n A m -i, n = 
2m ki A mn -\- (m — n — 3) A m -i, n+i 

+ (m + n + 3)A m+ i, n+1 + 2£ w ^ m ,„+i = 
2»M»+ (w + " — 3) A m -i, n -i 

+ (w — ;/ + 3)/l m+ i,„_i + 2k m Am,*-! = 

By means of these relations, all products can be expressed in terms of Aqq, 
Am, Aoi, and An- The following specific results are tabulated: 

A20 = i(4oo — 2kiAu — 2k Aio) I 
An = rr (kiAoo — 2A n — 2& ^oi) 

iki j 



(3.13) 



(3.14) 






(3.15) 



(3.16) 
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An = \{£u — k\Am — k>Au) 

An = -tj- (kiAio — Aoi — hAu) 

A30 = —£0^20 ~ kiAn 

Atx = — r (koAm + Au) 
ki 

The third-order product A21 is of considerable importance in the design 
of carrier amplifiers and radio transmitters, since the (2p — ^-product is 
the cross-product of lowest order falling back in the fundamental band when 
overload occurs. Figure 14 shows curves of A 2 i calculated by Mr. J. O. 
Edson from Eq. (3.12) by mechanical integration. 

We point out also that the linear-rectifier coefficients give the Fourier 
series expansion of the admittance of a biased square-law rectifier when two 
frequencies are applied. 

We shall next discuss the problem of reduction of the integrals appearing 
above to a closed form in terms of tabulated elliptic integrals . This can 
be done for all the coefficients above except the d 7 c for the zero-power law 
and for the d-c and two fundamentals for the linear rectifier. These contain 
the integral 

E(£o , ki) = I arc cos (k + h cos y) dy (3.17) 

Je 

which has been calculated separately and plotted in Fig. 22. When the 
arc cos term is accompanied by cos my as a multiplier with w^0,an integra- 
tion by parts is sufficient to reduce the integrand to a rational function of 
cos y and the radical \/l — (h + ki cos y) 2 , which may be reduced at once 
to a recognizable elliptic integral by the substitution z = cos y. It is 
found that all the integrals except that of (3.17) appearing in the results 
can be expressed as the sum of a finite number of integrals of the form: 

/cos 9 2 m j z 

By differentiating the expression z m 6 \/(l — z) 2 [l — (&o + h z 2 ] with 
respect to s, we may derive the recurrence formula : 

1_ 
(m - l)k\ 

+ (m- 2){kl - k\ - l)Z m _ 2 (3.19) 

— (2m — 5)k k 1 Z m - 3 + (m — 3)(1 — £o)Z m _ 4 ] 

9 Power series expansions of coefficients such as treated here have been given by A. G. 
Tynan, Modulation Products in a Power Law Modulator, Proc. I. R. E., Vol. 21, pp. 
1203-1209, Aug. 1933. 
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It thus is found that the value of Z m for all values of m greater than 2 can be 
expressed in terms of Z , Z h and Z 2 . 
Eq. (3.18) may be written in the form: 



if 



z m dz 



The substitution 



■\/(z — Zi) (z — z 2 )(z 3 — z)(z4 — z) 
zi= - (1 + k )/ki , % = - 1 

((1 - k Q )/ki , Case I) \ 
1, Case II / 

/ 1, Case I \ 

\ (1 - k )/k u Case II, / 

z 2 (z 3 — zi) — zi(z 3 — z 2 )u 



2 4 = 



z = 



z 3 — Zi — (z 3 — z 2 )« 2 



reduces the integral to 



^m 



*lV(Z4 - Z 2 )( 



Z3 



Zl) ° V(l - « 2 )(1 - ^W 2 ) 



where: 



V = 



Z3 — z 2 
Z3 — Zi 



2 _ (z 4 — zi)(z 3 — z 2 ) 

3/ ^— 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(Z 4 — Z 2 )(z 3 — Zi) 

Hence if K, E and II represent respectively complete elliptic integrals of 
the first, second, and third kinds with modulus k, and in the case of third 
kind with parameter — r\, we have immediately: 

2K 



Z = 
Zi = 



z, = 



klV(Z4 - Z 2 )(z 3 - Zl) 

2[z x K + (z 2 - z t ) n] 
hy/{zi — z 2 )(z 3 — Zi) 

-— / ==J=== \z\ K + 2z 1 (z 2 - z,)n 

ki V (Z4 — Z2) (Z3 — Zi) L 

2 r 1 rf« 1 

+ ^ - «o J (i _ , M 2 )V(i-« 2 )(i-^«')J 



(3.26) 
(3.27) 

(3.28) 
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To complete the evaluation of Z 2 , assume a relation of the following type 
with undetermined constants Ci, C?, G, C4: 

r du r du 

Jo (1 - nil 2 ? V(l - « 2 ) (1 - k 2 « 2 ) "" ! A V(l - « 2 ) (1 - k 1 « 2 ) 



V(i - « 2 ) 1 



Z\/(l - Z 2 )(l - K 2 Z 2 ) 



+ C 4 



1 - T/Z 2 



(3.29) 



u 

z 1.2 

u • 
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RATIO OF BIAS TO LARGER FUNDAMENTAL 



Fig. 15. — Fundamentals and (2p ± q) — product from full-wave biased zero-power-law 
rectifier with ratio of applied fundamental amplitudes equal to 0.5. Fi = larger funda- 
mental, F2 = smaller fundamental, F 3 = (2/> ± q) — product. 

Differentiate both sides with respect to z, set z = 1, and clear fractions. 
Equating coefficients of like powers of z separately then gives four simul- 
taneous equations in C\, d, Ca, C 4 . Solving for Cj, C 2 , C 3 and setting z = 1 
in (3.29) gives 



I (1 - 



du 



U 



K + 



ft (1 ~ T?« 2 ) 2 V(1 - «') (1 - K 2 « 2 ) 2(V 

(2 V - 3) k 2 - njr, - 2) "I 



2 



+ 



(3.30) 









164 



BELL SYSTEM TECHNICAL JOURNAL 



hi 

o _i 

? <■ 



i\ 
M 

o a. 
a. u 
q. o 

a: 



1.0 



0.8 



== Z 0.6 



0.4 



0.2 



-0.2 






Q2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

k 
RATIO OF BIAS TO L* a ^ c ° FUNnAMFW"*' 

Fig. 16. — Fundamentals and (2p ± q) — product from full-wave biased zero-power-law 
rectifier with equal applied fundamental amplitudes. 
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Fig. 17— The integral Z m with ki = 0.5. 

Since the necessary tables of II are not available, we make use of Legendre's 
Transformation, 10 which in this case gives: 



10 Legendre, Traitds des Fonctions Elliptiques, Paris, 1825-28, Vol. I, Ch. XXIII. 
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Fig. 18. — D-c. term in linear rectifier output with two applied frequencies. 
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Fig. 19. — Smaller fundamental in biased linear rectifier output. 



n = K + 



tan <}> 



775 t . ttk ^E{4>) - BF(4>)] 
yl-r sim <p 



1/2 

<f> = arc sin — 

K 

dd 



F{<t>) = [ r/r -> n 

Jo yl -k sin 2 6 

r* 

E(<t>) = / Vl - k 2 sin 2 d6 
Jo 



(3.31) 
(3.32) 
(3.33) 
(3.34) 



The functions F(#) and E(<£) are incomplete elliptic integrals of the first 
and second kinds. They are tabulated in a number of places. Fairly good 
tables, e.g. the original ones of Legendre, are needed here since the difference 
between KE(0) and EF(<£) is relatively small. 
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Summarizing: 



Fig. 22.— Graph of the integral E (*o. ki). 



Case I, * + *i > 1, *o- *i < 1 

Zo ~ V*7 

Zi = ^ [££(*) - EF(4>)] - Z ( , 



il 1 



Z 2 = t I z - hZ x - 



Vki 



-A 



y ■ 4*1 



= arc sin 



2*i 



+ *o + *i 

Case II, h + *i < 1, *o - *i > -1 
2tf 



Zn = 



V(i + *0 2 - * 



Z x = 4 lKE(<f>) - EF(<f>)} - Z„ 



(3.35) 



Zt = Tp |(1 + k\ - kl)Z - 2k kiZi - IE V(l + *i) 2 - ttj 



■ (3.36) 



4*, 



<f> = arc sin 



(1 + hY - *5 



— k + *i 
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The values of the fundamentals and third-order sum and difference 
products for the biased zero-power-law rectifier have been calculated by the 
formulas above for the cases k\ = .5 and k\ = 1. The resulting curves are 
shown in Fig. (15) and (16). The values of the auxiliary integrals Zo, Z\ , 
and 7,i are shown for k\ = .5 in Fig. (17). These integrals become infinite 
at ko = 1 — ki so that the formulas for the modulation coefficients become 
indeterminate at this point. The limiting values can be evaluated from 
the integrals (3.3), etc., directly in terms of elementary functions when the 
relation k = 1 — k\ is substituted, except for the E-function. 

Limiting forms of the coefficients when k is small are of value in calcu- 
lating the effect of a small signal superimposed on the two sinusoidal com- 
ponents in an unbiased rectifier. By straightforward power-series expan- 
sion in ko, we find : 
Zero-Power-Law Rectifier, k Small: 



A 1Q = —E - 



2E 



7T 2 (1 - k 2 ) 



kl + 



A 01 = -i_ [ E - (1 - k\)K] + 2 



^21 = 



7T 2 £, 

4 
37T 2 £l 



E 



7T 2 k\ \1 — k\ 

[(1 - 2k\)E - (1 - k\)K] 



-K)k- + 



+ ^[*-^H*» + 



(3.37) 



In the above expressions, the modulus of K and E is k\. When k = 0, 
these coefficients reduce to half the values of the full-wave unbiased zero- 
power-law coefficients, which have been tabulated in a previous publication. l 
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